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Abstract: 

We present a calculation of next-to-leading order QCD corrections to jjj production via 
weak-boson fusion at a hadron collider in the form of a flexible parton-level Monte Carlo program 
which allows us to study cross sections and distributions within experimentally relevant selection 
cuts. The radiative corrections are found to be moderate with residual scale uncertainties being 
considerably improved beyond tree level. The sensitivity of the reaction to anomalous W + W~j 
couplings is investigated. 



1 Introduction 



After the start-up of the CERN Large Hadron Collider (LHC) weak-boson fusion (WBF) reac- 
tions deserve attention as a powerful means for discovering the Higgs boson predicted by the 
Standard Model (SM) and for a later determination of its properties [lr[7]. Moreover, WBF 
processes will shed light on the mechanism of electroweak symmetry breaking itself, revealing 
whether the SM is truly realized in nature or extensions are necessary. No matter which scenario 
proves suitable for describing the physics of elementary particles, WBF reactions are expected to 
provide valuable information on the interactions among weak gauge bosons. Because of the strik- 
ing kinematic features of this class of processes with well separated tagging jets in the forward 
and backward regions of the detector, a priori overwhelming QCD backgrounds can be removed 
efficiently from the electroweak signal signature [3,4,8]. Radiative corrections have been shown 
to be well under control [9l-ll4j. and the interference of the purely electroweak WBF processes 
with QCD production modes giving rise to the same final-state configurations has been found 
to be negligible [JTJ[T5l[16] ■ Various WBF processes have thus also been considered as possible 
means for accessing couplings of the Higgs boson to weak gauge bosons [6], triple and quartic 
weak boson vertices [H1[T8], both in the context of the SM and possible extensions thereof. 

In this work, we consider the electroweak (EW) production mode pp — > jjj at order C(a 3 ), 
which provides access to the coupling of the photon to weak gauge bosons and may thus yield 
information on this triple gauge boson vertex complementary to bounds derived mainly from 
single- and double gauge boson production at LEP and the Tevatron [19]. We develop a fully 
flexible parton-level Monte Carlo program, structured similarly to an existing code |20| for re- 
lated WBF reactions. The program allows for the calculation of cross sections and distributions 
within experimentally feasible selection cuts. Interference effects that are entirely negligible 
in the phase-space regions where WBF can be observed experimentally [l H ll5 t fT6] are disre- 
garded throughout. In this approximation, the next-to-leading order (NLO) QCD corrections 
to pp — > jjj are computed. In order to obtain well-defined predictions for the photon+2 jet 
final state, we employ the isolation criterion of Frixione |21j . In addition to discussing the SM 
production process and estimating the theoretical stability of the prediction, we investigate the 
sensitivity of the reaction to the W + W~ r y vertex. Although our discussion focusses on jjj pro- 
duction at the LHC, our results can readily be applied to any other high-energy hadron-hadron 
collider such as the Tevatron. 

We start with a brief outline of the calculational techniques used for developing a parton-level 
Monte Carlo program in the context of the SM and a summary of the consistency checks we have 
performed on the thus obtained code in Sec. [2j Numerical results are presented and discussed in 
Sec. El Our conclusions are given in Sec. [H 

2 Framework of the calculation 

At tree-level, EW jjj production in hadronic collisions mainly proceeds via quark scattering, 
qq' — > jqq', mediated by the exchange of a weak gauge boson or a photon. The relevant charged- 
current (CC) Feynman-diagrams can be classified in terms of the two topologies depicted in Fig. [1] 
for a specific subprocess, depending on whether the photon is being emitted off the t-channel 
exchange boson (a) or off a fermion line (b). To neutral-current (NC) jjj production, only graphs 
of the latter topology contribute. Diagrams for the related anti-quark initiated subprocesses 
qq' — > qq'j and qq' — > qq'j are easily obtained thereof via crossing. The evaluation of the 
relevant tree-level contributions is performed by means of the amplitude techniques of Ref. |22j . 
In the phase-space regions where WBF can be observed experimentally, with widely separated 
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Figure 1: Feynman-diagram topologies contributing to the Born process us — > dcj. Graphs 
analogous to (b), with the photon being emitted off the lower quark line, are not shown. 



jets of large invariant mass, contributions from annihilation processes such as qq' — > W* — > W7, 
with subsequent decay W — > qq' , are negligible. Identical- flavor interference effects, which require 
the interference of i-channel with u-channel diagrams, are strongly suppressed both at tree-level 
and beyond in WBF-type reactions, as explicit calculations for the Higgs-production mode have 
revealed [11|J16|. We therefore entirely disregard these types of contributions. In the following, 
we will refer to pp — > jjj within the mentioned approximations as "EW jjj production" . 

The structure of the NLO QCD corrections to EW jjj production is very similar as in the 
case of related WBF processes. We therefore do not present details of the calculation reminiscent 
of what has already been described in the literature [9lll2p i3j. but focus on aspects specific to the 
jjj final state. The real-emission contributions to the CC and NC production modes discussed 
above are obtained by attaching a gluon to the (anti-)quark lines in all possible ways. Crossing 
these diagrams yields contributions with an initial-state gluon and an extra anti-quark in the 
final state, such as gq' — > jqqq' ■ Singularities in the soft and collinear regions of phase space 
are regularized in the dimensional-reduction scheme with space-time dimension d = 4 — 2e. 
The cancellation of these divergences with the respective poles of the virtual contributions is 
performed by introducing the appropriate counter-terms of the dipole-subtraction method [23] , 
Since the color structure of the reaction under consideration is identical to the related EW Zjj 
process, the respective counter terms can be straightforwardly adapted from Ref. [12j. 

The virtual contributions are obtained from the interference of the sum of all one-loop di- 
agrams, .My, with the Born amplitude, Mb- Within the approximations outlined above only 
selfenergy, vertex, and box corrections to either the upper or the lower quark line have to be con- 
sidered, while contributions from diagrams where a gluon is exchanged between the two fermion 
lines vanish, due to color conservation. In order to compute the relevant diagrams, we have 
separately calculated all one-loop corrections to a quark line with only one t-channel exchange 
boson (V = W^, Z,j) being attached, 

q(ki)^V( gi )q(k 2 ), (1) 

and to a quark line with the t-channel exchange boson and the external photon being attached, 

q(k l )^V(q l ) 1 (q 2 )q(k2). (2) 
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The calculation of the virtual contributions is performed analytically by means of Mathematica, 
making use of the tracer package |24| . 

The interference of all diagrams containing one-loop corrections to either the upper or the 

(i) 

lower quark line, Ai v , with the entire Born amplitude is of the form 
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where Qi is related to the momentum transfer between the respective initial- and the final-state 
quarks via Q\ = —{k\ — /C2) 2 , is the renormalization scale, Cf = 4/3, c v i rt = 7r 2 /3 — 7, and 



Aiy is a finite remainder. 

The poles in Eq. (|3|) are canceled by respective singularities in the phase-space integrated 
counter terms, which in the notation of Ref. |23| are given by 
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The non-factorizable terms Aiy can be expressed in terms of the finite parts of the Passarino- 
Veltman Bq, Cjk, and Djk functions, which are evaluated numerically. 

In order to ensure the reliability of our calculation, several checks have been performed: 

• The tree-level and real-emission amplitudes have been compared to the corresponding ex- 
pressions generated automatically by MadGraph [25j for a representative set of phase-space 
points. We found full agreement. 

• The integrated LO cross sections have been compared to the corresponding results of the 
MadEvent package [26] for inclusive selection cuts. The cross sections agree within the 
numerical accuracy of the two programs. 

• To validate our implementation of the dipole subtraction, we checked that the real-emission 
contributions approach the subtraction terms in singular regions of phase space. 

• The real-emission matrix elements have been found to vanish when the polarization vector 
of the external gluon is replaced with its momentum. This procedure tests the QCD gauge 
invariance of the real-emission amplitudes. 

• The tree-level, real-emission, and virtual matrix elements have been found to vanish when 
the polarization vector of the external photon is replaced with its momentum. This proce- 
dure tests the QED gauge invariance of the amplitudes. 

• In order to test our virtual corrections at an analytical level, we have expanded our 
Mathematica modules for the q — > V^q building blocks of Eq. ([2]) to allow for the si- 
multaneous evaluation of q — > VZq. The latter process is related to qq — > VZ via crossing. 
Comparing our thus obtained virtual corrections for the qq — > ZZ mode to the litera- 
ture |27j at an analytical level provides a valuable check of our calculation. 

• In complete analogy to the jjj case, we have determined the real and virtual NLO-QCD 
corrections for EW Zjj production, implemented them in our Monte-Carlo program and 
compared to the code of Ref. |12j . Again, we found full agreement within the numerical 
accuracy of the two codes, both at amplitude level for selected phase-space points and for 
integrated cross sections. 
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3 Numerical results and discussion 



The cross-section contributions discussed above have been implemented in a fully flexible parton- 
level Monte-Carlo program which allows us to compute cross sections and distributions within 
experimentally relevant selection cuts. 

For our numerical analysis we employ the parton distribution functions of the CTEQ collab- 
oration with parameterizations being provided both at LO and NLO QCD. Specifically, we use 
the CTEQ6M set with a s (m z ) = 0.118 at NLO, and the CTEQ6L parameterization at LO [28] . 
For reference, we have also calculated cross sections with the MSTW parton distributions |29j . 
which are available only at NLO QCD. We found that the respective results differ from those 
obtained with the CTEQ6 distributions by less than 2%. Since in our calculation quark masses 
are neglected, we entirely disregard contributions from external b and t quarks. As electroweak 
input parameters, we choose the weak boson masses, mz = 91.188 GeV, mjy = 80.423 GeV, and 
the Fermi constant, Gf = 1.166 x 10~ 5 /GeV 2 . The other parameters, oqeDi and sin 2 #w, are 
computed thereof via tree-level electroweak relations. Final-state partons are recombined into 
jets according to the hr algorithm [30] with a resolution parameter D = 0.7. If not indicated 
otherwise, results are shown for pp collisions at the LHC design energy of v 7 ^ = 14 TeV. For the 
Cabibbo-Kobayashi-Maskawa matrix, VckMj we have used a diagonal form, equal to the identity 
matrix, which is equivalent to employing the exact VckMj when the summation over final-state 
quark flavors is performed and when quark masses are neglected, cf. Ref |14| . 

For our phenomenological analysis we utilize selection cuts motivated by Ref. [I] to allow for 
the identification of WBF events with photons in the final state at the LHC. We require at least 
two hard jets with 

p Tj >20GeV, y<4.5, (5) 

where pTj denotes the transverse momentum and yj the rapidity of a jet j being reconstructed 
from massless partons of pseudo-rapidity \r]j\ < 5. The two jets of highest transverse momentum 
are referred to as "tagging jets" . We impose a large rapidity separation between the two tagging 
jets, 

A vjj = Ivji -vd > 4 - 4 > (6) 

and furthermore demand that they be located in opposite hemispheres of the detector, 

Vji x Vj2 < , (7) 

with an invariant mass 

Mjj > 600 GeV . (8) 

Imposing the jet-selection criteria of Eqs. and additional cuts on the final-state photon 

to be discussed below, the LO differential cross section for EW jjj production is finite. At 
NLO, initial-state singularities due to collinear q — ^ qg and g — > qq splittings can arise. These 
are taken care of by factorizing them into the respective quark and gluon distribution functions 
of the proton. Additional divergences stemming from the t-channel exchange of low-virtuality 
photons in real-emission diagrams are avoided by imposing a cut on the virtuality of the photon, 
min = 4 GeV 2 . Events that do not satisfy the Q 2 min constraint would give rise to a q — > q*y 
collinear singularity, which is part of the QCD corrections to p7 — > jjj and not taken into 
account here. We have checked that the NLO-QCD cross section within typical WBF cuts is 
quite insensitive to the virtuality cutoff, changing by less than a permille when Q 2 min is lowered 
from 4 GeV 2 to 0.1 GeV 2 . 
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To ensure that the final-state photon is well-observable in the central-rapidity range, we 
require its transverse momentum pry and rapidity y 7 to fulfill the following conditions: 

PT 7 >20GeV, |y 7 |<2.5, (9) 

min(?/ji,?/j2) +0.7 < y 7 < max(y jl ,y j2 ) - 0.7. (10) 

In order to isolate the photon in an infrared-safe way from partons without having to introduce 
parton-to-photon fragmentation contributions, we apply the criterion suggested in Ref. [21j : An 
event is considered as acceptable only, if the hadronic energy deposited in a cone around the 
direction of the photon is limited by 

£ PTi < !~ cos V 7 (vi?<5 ). (ii) 

Here, the summation index i runs over all final-state partons found in a cone of size R in the 
rapidity-azimuthal angle plane around the photon, pTi denotes the transverse momentum, and 
i?i 7 the separation of parton i from the photon, while Sq stands for a fixed separation. To 
illustrate the isolation-cone dependence of the integrated cross section for EW jjj production 
within the cuts of Eqs. (f5])~ (|lip . we list <r cuts for various values of 5q in Table [TJ At LO with 
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Table 1: EW pp — > jjj cross sections at the LHC with y/S = 14 TeV at LO and NLO within the 
cuts of Eqs. (j5j)~ (fTTj) for fip = [ir = Qi and four different values of <5o- The relative statistical 
errors of the quoted results are at the sub-permille level. 

only two partons in the final state which, due to the selection cuts we impose, are typically 
well-separated from the photon, the 5q dependence is very small, amounting to less than 4% 
in the considered range. A slightly larger isolation-cone dependence occurs in the NLO cross 
sections, as the additional final-state parton in the real-emission contributions can come close to 
the photon and therefore is more sensitive to the isolation criterion. Alternatively to the photon 
isolation procedure of Ref. [21], we could opt for computing the full photonic final state, including 
contributions which stem from parton-to-photon fragmentation. This procedure would allow us 
to drop the criterion of Eq. (|lip . but at the same time introduce a dependence on the so far 
only poorly known photon fragmentation functions [H]. In the following we therefore restrict 
ourselves to the direct photon production contributions, removing all collinear fragmentation 
contributions with the help of the isolation procedure described above, and set 5q equal to 1.0. 

To illustrate the impact of Eq. (jlip on <j cuts , in Fig. [2] we show the rapidity-azimuthal angle 
separation of the photon from the closest jet, da/dRJ^ n , at LO and NLO. Events with a photon- 
jet separation smaller than 5o yield non-vanishing contributions, unlike what would be observed 
in the case of an explicit cut on the observable, such as i?j 7 > 5o- We furthermore observe that, 
in complete analogy to the case of EW Zjj production discussed in Ref. |12] , the real-emission 
contributions shift da/dR™ n to smaller values, thereby reducing the isolation of the photon from 
the jets. 

In order to assess the dependence of <j cuts on unphysical scales, in Fig. [3] we display the 
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Figure 2: Rapidity- azimuthal angle separation of the photon from the closest jet in EW pp — > jjj 
production at the LHC with ^/S = 14 TeV at LO (dashed black line) and NLO (solid red line) . 
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Figure 3: Dependence of the EW pp — > -yjj cross section at the LHC with \/~S = 14 TeV on the 
factorization and renormalization scales. The NLO curves show cr cuts as a function of the scale 
parameter ^ for three different cases: fiR = fiF = £,Qi (solid red line), [If = £Qi and fiR = Q{ 
(dot-dashed blue line), [Ifi = £Qi and fiF = Qi (dashed green line). The LO cross sections 
depend only on fip (dotted black line) . 
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Figure 4: Transverse momentum distribution of the highest-p-p tagging jet in EW pp — > jjj 
production at the LHC with \fS = 14 TeV at LO (dashed black line) and NLO (solid red line) 
[panel (a)] and relative corrections according to Eq. ()13f) when the factorization and renormal- 
ization scales are varied in the range Qi/2 < /zr = ^f < 2Qi [panel (b)]. 



integrated cross section as a function of the renormalization and factorization scales, /Ur and hf, 
which are taken as multiples of the scale parameter /j,q, 

Mf = £f Mo , MR = £r Mo , (12) 

where for each fermion line /io is identified with the respective momentum transfer Qi between 
the incoming and the outgoing parton. The LO cross section only depends on fip. For the central 
scale, = fJ'R = Qi, the NLO-QCD corrections amount to 12%, being slightly larger than in 
the case of massive gauge-boson production in WBF-type reactions [12]. Qualitatively similar 
results are obtained for alternative choices of the factorization and renormalization scales such 
as the average transverse momentum of the iV jets in an event, ^PTj/N . In the following, we 
set fiF = Mr = Qii unless stated otherwise. 

In addition to the normalization of the integrated cross section, radiative corrections can 
affect the shape of various kinematic distributions. To quantify the relative size of NLO-QCD 
corrections to an observable together with the residual scale uncertainties, we introduce the 
quantity 5(0), defined by 

x(fn\ da(nF,HR)/dO 

6{0) = d *^o {flF = ^ = Qt)/d0 -^ 

as a measure for the deviation of the LO or NLO expression da(np, Mr), evaluated at an arbitrary 
scale, from the corresponding NLO observable at the default scale fXp = /iR = Qi. 

Figured] shows the transverse momentum distribution of the highest-pr tagging jet at LO and 
NLO for our default settings, together with the associated relative corrections for two different 
values of the scale parameter £ = £r = £f = 1/2 and 2. The difference between the curves 
for the two values of £ indicates the scale uncertainty of da/dpxj at LO (dashed black lines) 
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Figure 5: Transverse momentum distribution of the highest-p-p tagging jet in EW pp — > jjj 
production at the LHC with \/S = 7 TeV at LO (dashed black line) and NLO (solid red line) 
[panel (a)] and relative corrections according to Eq. ()13f) when the factorization and renormal- 
ization scales are varied in the range Qi/2 < ^ = < 2Q, L [panel (b)]. 



and NLO (solid red lines), respectively. At low transverse momenta, relatively large positive 
corrections occur. Towards higher values of pTj the impact of NLO-QCD corrections on the 
cross section decreases. In this kinematic range the scale uncertainty of the LO distribution 
is sizeable, but can be efficiently reduced by taking the NLO-QCD contributions into account. 
Similar features are observed for lower values of the hadronic center-of-mass energy, cf. Fig. [5J 
Decreasing from 14 TeV to 7 TeV yields smaller cross sections and slightly larger scale 
uncertainties with radiative corrections again affecting predominantly the low-p^- region. 

Numerically, the EW jjj cross section is dominated by CC subprocesses, where the two 
incoming partons scatter by t-channel VF-boson exchange. It is thus particularly suitable to 
access the W + W~^ coupling. Since NLO-QCD corrections do not affect the weak boson vertex, 
the implementation of anomalous triple gauge boson couplings (TGC) in our Monte-Carlo pro- 
gram is straightforward, allowing us to study their impact on the cross section and kinematic 
distributions. In the effective Lagrangian approach of Ref. [32], the Lorentz structure of the 
momentum-space vertex W~ (q)Wt (q')'j l j i (p) (with all momenta outgoing, p + q + q' = 0) is 
parameterized in the form 



^Jq, q',p) = q' a g^ ( 2 + Are 7 + A 7 — ^-^ - gV + Are 7 + A 7 ^) 
+ {q'^-q^) 



2 m l w ) m A w 



(14) 



where C and P invariance have been imposed and only operators up to dimension six are taken 
into account. For an on-shell photon we have p 2 = 0, and thus the third term in the above 
expression further simplifies. In the limit Are 7 — > 0, A 7 — > 0, Eq. (]14p reduces to the SM 
expression for the W + W~^f vertex, which is compatible with combined limits from LEP and the 
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Figure 6: Transverse momentum distribution of the photon in EW pp —> production at the 
LHC with \/S = 14 TeV at LO (dashed black line) and NLO (solid red line) [panel (a)] and 
relative corrections according to Eq. (|13|) when the factorization and renormalization scales are 
varied in the range Qi/2 < pr = pp < 2Qi [panel (b)]. 



Tevatron [33J. A priori, the effective Lagrangian approach gives rise to unitarity violations in the 
high-energy regime. To avoid this unphysical behavior, the anomalous gauge boson couplings in 
Eq. (fT4j) have to be supplemented by a form factor. Following the prescription of Ref. [17] . we 
perform the substitutions 



Ak 1 



Ak 7 



A 7 



A" 



1 + 



A 5 " 



1 + 



1 + 



l<? 2 l 



(1 + 



A- 
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where n is an exponent that has to be chosen large enough to ensure unitarity and A is interpreted 
as the scale where new physics enters that is not accounted for by the effective Lagrangian 
approach. Below, we will set A = 2 TeV and n = 1. We will show, however, that the impact of 
the form factor on the observables we consider is rather small. 

Anomalous gauge boson couplings mainly affect events at high energy or large transverse 
momenta. Therefore, the high-p^-y tail of the photon's transverse momentum distribution is 
expected to be particularly sensitive to Ak 7 and A 7 . For being able to distinguish the anomalous 
coupling effects from radiative corrections, a precise knowledge of this distribution both within 
the SM and its considered extension is essential. The SM prediction, da/dpx-y, is displayed 
at LO and NLO together with the associated relative corrections for different values of the 
factorization and renormalization scales in Fig. [6j The NLO-QCD corrections to this distribution 
are moderate, but distort its shape noticeably. Reminiscent of our observations for the transverse 
momentum distribution of the high-py tagging jet, we find that in the low transverse momentum 
range relatively large positive corrections occur. The scale dependence increases with px-y, but 
is considerably improved when going from LO to NLO. Having estimated the uncertainties of 
the SM distribution, we now turn to an analysis of the impact anomalous TGC may have on 
da I dpTy ■ 
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Figure 7: Ratio of da^-P /dpr^ and da^° /dpr^ in EW pp — > jjj production at the LHC with 
y/S = 14 TeV for Ak 7 = 0.02, A 7 = (dashed green lines), for Ak 7 = -0.02, A 7 = (solid green 
lines), for Ak 7 = 0, A 7 = 0.02 (dashed blue lines), and for Ak 7 = 0, A 7 = -0.02 (solid blue 
lines). In panel (a) the form factor of Eq. (|15p is included, while in panel (b) no form factor has 
been supplemented. 



To this end, we consider the ratio of the distribution with non- vanishing anomalous TGC, 
da^}P / dpT-y , to the respective SM curve, da^° /dpr-y, depicted in Fig. [7] for two scenarios: In 
panel (a) the form factor of Eq. (|15p is included, while in panel (b) no form factor has been 
supplemented. The two parameters Ak 7 and A 7 are set to non-zero values compatible with 
current experimental limits. We find that da/dpx^ is sensitive to A 7 over the entire transverse 
momentum range. As expected, a non-zero value of A 7 becomes particularly evident towards 
the high-pr 7 tail. The dependence on Ak 7 is less pronounced. The impact of the form factor 
is small, as becomes evident from contrasting the results in panel (a) with the respective curves 
in panel (b). Only at very large values of pT-y the distributions without a form factor start to 
exceed the corresponding predictions including a form factor. We note that anomalous TGC give 
rise to similar shape distortions in the tagging jets' transverse momentum distributions. 

Further information on the structure of the W + W~~f vertex can be obtained from the 
azimuthal-angle separation of the two tagging jets, Acftjj. Azimuthal-angle correlations have 
been suggested as a tool for determining the CP properties of the Higgs boson in pp —> Hjj [5]. 
Eboli and Garcia [17] have furthermore considered da/dA(j)jj as a means for constraining anoma- 
lous TGC in pp — > W jj. Figure [8] shows this distribution, normalized to the SM result, for 
various values of Ak 7 and A 7 in pp — > jjj. A marked enhancement is observed at low values of 
A(f)jj for non-zero values of A 7 , while the impact of Ak 7 is very small. Including the form factor 
of Eq. (fT5j) does not change the qualitative features of the distribution. 
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Figure 8: Ratio of da^° /dAfaj and da^° /dAfaj in EW pp — > 7jj production at the LHC 
with y/S = 14 TeV for = 0.02, A 7 = (dashed green lines), for Ak 7 = -0.02, A 7 = (solid 
green lines), for Ak 7 = 0, A 7 = 0.02 (dashed blue lines), and for Ak 7 = 0, A 7 = —0.02 (solid 
blue lines). In panel (a) the form factor of Eq. (|15p is included, while in panel (b) no form factor 
has been supplemented. 

4 Conclusions 

In this article we have presented results for EW jjj production in hadronic collisions at NLO- 
QCD accuracy, obtained with a fully-flexible parton-level Monte Carlo program which allows 
us to study cross sections and distributions within realistic experimental selection cuts. NLO- 
QCD corrections were found to enhance integrated cross sections by about 12% and change 
the shape of some distributions. The residual scale dependence of all results is significantly 
improved beyond the Born approximation. After the inclusion of NLO-QCD corrections, the 
cross section within the WBF-specific selection cuts we have introduced in Sec. [3] changes by 
less than 2% when factorization and renormalization scale are varied simultaneously in the range 
Qi/2 < [J>f = [J-R < 2Qj. Compared to the more pronounced scale dependence of about 9% at 
LO, this improvement indicates that the process is under excellent control perturbatively. 

Moreover, we found that pp — > jjj could yield information on the VF + VF~7 coupling com- 
plementary to bounds derived from single and double gauge boson production at LEP and the 
Tevatron. To demonstrate the sensitivity of the EW 7jj cross section on this vertex, we im- 
plemented anomalous TGC in our code and demonstrated that noticeable changes are induced 
by non-SM contributions to the M /+ W /_ 7 vertex on selected observables, such as the transverse 
momentum distributions of the photon and of the tagging jets, as well as the azimuthal angle 
separation of the tagging jets. Since anomalous gauge boson couplings affect the shape of these 
distribution in a manner very different from NLO-QCD contributions, their sensitivity is not 
spoiled by radiative corrections. 
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